An Explanation for Heavy Quark Energy Loss Puzzle by Flow Effects 
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Based on our new model potential in the presence of collective flow describing the interaction 
of the hard jet with scattering centers, an explanation for heavy quark energy loss puzzle is given 
by collective flow effects. It is shown that the collective flow changes the pole of propagator which 
lead to reduce significantly the dead cone and change the LPM destructive interference comparing 
to that in the static medium. Consider the collective flow with velocity v z along jet direction, the 
heavy quark energy loss is nearly the same as light quark in the region of 0.065 < v z < 0.13. It turns 
out that the average v z of rescattering centers is around 0.08, the differences of the effective average 
energy loss among charm, bottom and light quarks are very little from a full 3D ideal hydrodynamic 
simulation for 0-10% central events of Au-Au collision at RHIC energy by including collective effects. 
The flow dependence of the quark energy loss implies that the light and heavy quark have almost 
the same suppression of high pr hadron spectrum in high-energy heavy-ion collisions. 

PACS numbers: 12.38.Mh, 24.85.+p; 25.75.Bh; 25.75.-q 



Introduction — The formation and observation of 
quark-gluon plasma (QGP) are the most important goals 
in high energy heavy ion physics. Jet quenching [l| or 
suppression of large transverse momentum pt hadrons 
caused by the energy loss of a propagating parton in a 
dense medium, has become a powerful tool for the study 
of properties of QGP. Light quark and gluon jet quench- 
ing observed via %, rj suppression in Cu+Cu and Au+Au 
collisions at = 62 — 200GeV at the relativistic heavy 
ion collider (RHIC) has been remarkably consisted thus 
far with predictions @, d, 0, Q from jet quenching theory. 
It is believed that a good test of QGP formation is the 
open charm suppression, which can be measured at RHIC 
by comparing the pt distributions of single electron pro- 
duction in d+Au and Au+Au collisions. It is estimated 
that the heavy quark mass leads to a kinematical "dead 
cone" effect that reduces significantly the induced radia- 
tive energy loss of heavy quarks [E 0, H, HI based on the jet 
quenching picture. However, recent non-photonic single 
electron data [HI, [H, 13 shows that a much larger 
suppression of electrons, which is similar to the suppres- 
sion of light hadrons, than predicted was observed in the 
Pt = 4 — 8 GeV region. Later the collisional energy loss 
is included to explain the non-photonic single electron 
data, but this still doesn't yield a consistent description. 
Recently the ADS / CFT approach has been used for this 
question [l5j]. At present the heavy quark energy loss be- 
comes a puzzle to challenge the scenario of jet quenching. 

The currently available studies suffer from one cru- 
cial drawback, the medium induced radiative energy loss 
should be computed in a expanding QCD medium with 
collective flow instead of a static medium. The medium 
produced in nucleus-nucleus collisions at RHIC equili- 
brates efficiently and builds up a flow field. This flow 
changes the color-electric field and produces a color- 
magnetic field [j^ j. and will surely change the parton 
energy loss. In this letter, based on our previous new 
model potential [l6| with collective flow, we will report a 
first study of the heavy quark energy loss in the presence 



of collective flow in perturbative Quantum Chromody- 
namics (pQCD) and give an explanation for heavy quark 
energy loss puzzle in the framework of jet quenching the- 
ory. Based on our new potential model with collective 
flow, we find that the flow effects reduces significantly 
the dead cone of heavy quark jets. We are led to the 
conclusion that the decrease of the dead cone increases 
the heavy quark energy loss. Using the full 3D ideal 
hydrodynamic simulations [l7j . we obtained the average 
flow velocity of rescattering centers and effective aver- 
age energy loss of heavy quark jets. It turns out that 
the effective average energy loss of heavy quark jets is 
nearly the same as the light quark jets, this agrees with 
the experimental data that the suppression of electrons 
for heavy quarks is similar to light hadrons. 

Flow Effect on Gluon Radiation for Heavy Quark 
Jet — In our previous work we model the jet 
interactions in the QGP in the presence of collective 
flow as random color screened potentials A"{qi) = 
{V/ low ( qi ),A{ l ™( qi )), 

V/ low ( qi ) = 27r%?-vq) e -^ x ^(q)T a (i)T a (i) , (1) 
Af™(%) = 27r ( 5(gO-v.q)ve-^ x ^(q)T a (j)T a (z),(2) 

where v(q) = 47ra s / (q 2 - (v • q) 2 + /i 2 ), q l = (<?", q) is 
the momentum transfer, fi the Debye screening mass, v 
the flow velocity, T a (j) and T a (i) the color matrices for 
the jet and target parton at position x», respectively. The 
collective flow changes the color-electric field, produces a 
color-magnetic field and leads to non-zero energy transfer 

= v ■ q from the target parton to the jet which differs 
with Gyulassy- Wang's static potential model 

Here we investigate the rescattcring-induccd gluon ra- 
diation of heavy quark jet by considering the flow effect 
resulting from the moving parton target. We will work in 
the framework of opacity expansion developed by Gyu- 
lassy, Levai and Vitev (GLV)[l|| and Wiedemann 
The opacity is defined as the mean number of collisions 
in the medium. In the framework of the opacity cxpan- 
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sion, the first order opacity contribution is dominant be- 
cause the higher order corrections contribute little to the 
radiative energy loss as shown by GLV [T^ |. 

In order to compute the medium-induced radiative en- 
ergy loss to the first order in the opacity expansion, 
we need to obtain the amplitudes for self-quenching in 
the vacuum, and for single rescattcring and double Born 
scattering with the same target in the QCD medium in 
the presence of collective flow. Consider a hard heavy 
quark jet with mass M and initial energy E produced 
at x = (z ,x ^). It interacts with the target parton 
at xi = (zi,xi±) with flow velocity v by exchanging 
gluon with four-momentum q, radiates a gluon with four- 
momentum k and polarization e(fc), and emerges with 
final four- momentum p. In Ref. [6j , it is shown that glu- 
ons in the medium can be be approximated as massive 
transverse plasmons with mass m g m fi/y/2. Then p, 
k and polarization e(fc) can be written in the light-cone 
components, 
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where u = xE, E + = 2E > p.. 

At zeroth order in opacity, the jet has no interaction 
with the target parton, we obtain the same radiation am- 
plitude as that in the static medium in Ref. 
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which indicates the depletion of radiation in the "dead 
cone" [1[ at angles 



8 <8 cl = Jm 2 +x 2 M 2 /(xE) . 



(7) 



Based on our new potential with the collective flow of 
the quark-gluon medium in Eqs. fl]) and ([2]), assuming the 
flow velocity |v| <C 1, the radiation amplitudes associated 
with a single scattering and the double Born rescattering 
with two scattering centers in the contact limit can be 
expressed as 

R<-^ = 2ig(c 2 {e l ^T Zw -e l "° + ^ M ^)[T C , T a ] 

+H 1 e i "° + ""r" M zio r e r a + Hi(l-« a ) 2 (e- i ^fr* 10 
-e r ° + "-^ M z ">)T a T c ) ■ 6 X (1 - v z ) , (8) 
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where z±o — z\ — Zo, Cr and Ca are the Casimirs of jet 
in fundamental representation in cIr dimension and the 
target parton in adjoint representation in cLa dimension, 
respectively. 
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Different from the static medium case, the single and 
double Born scattering amplitudes depends on the col- 
lective flow of the medium. Our results agree with the 
GLV results @ in static medium at zero flow velocity. 

In calculating radiation amplitude in Eqs.((8|) and (J9j) , 
compared to the static medium case we can see that the 
collective flow changes the poles of the propagator as 
follows: the poles of the quark propagators are shifted 
by a flow velocity-dependent term lom and then times a 
flow factor 1/(1 — v z ). The pole shift leads to that the 
mass-dependent energy difference oj m in Ref. @ shifts 
to ui m — ojm as shown in Eqs. ([5]) and ©. For the pro- 
cesses that heavy quark jet interacts directly with the 
target partons, the contribution of shifted poles of both 
quark propagators cancels each other, the net contribu- 
tion keeps the same radiation amplitude as that in the 
static medium except the phase factor, as shown by Hi 
and Ci terms in Eqs. JSJ) and ([9]). From Hi and Ci in 
Eqs. (TT2")) and (TTJ} we can see easily that these processes 
have the dead cone Q c \ in Eq. which is same as that 
in the static medium. For the processes that heavy quark 
jet interacts with target partons via triplc-gluon vertex, 
the contribution of shifted pole of the quark propaga- 
tor can't cancel the contribution coming from the pole 
of gluon propagator, the flow effects changes the radia- 
tion amplitude as shown by H 2 and C 2 terms in Eqs. 
© and ©. From H 2 and C 2 in Eqs. flT3J) and (HD we 
can see clearly that in these processes the collective flow 
decreases dead cone 9 c i as 



m 2 + x 2 M 2 - (2v z xM 2 )/(l - v z )/(xE) . (16) 



Compared to the dead cone term M 2 /E 2 coming from 
the mass effects, the dead cone change 2v z M 2 /[(l — 
v z )xE 2 ] coming from the flow effects is proportional to 
l/x = E/uj, this term deduces significantly the mass- 
effect dead cone for high energy heavy quark jet. The 
real dead cone with collective flow effects should be 



mm[e cl ,6 c2 ] = 9 c2 > 0. 



(17) 
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FIG. 1: The ratio of dead cone in the medium with and with- 
out flow as function of flow velocity in the jet direction when 
lo/E = 0.2 

Fig[T] is the ratio between the dead cone in the medium 
with flow and that in the static medium as function of 
flow velocity v z as ui/E = 0.2. It shows that the dead 
cone decrease to zero rapidly with increasing v z . The 
dead cone of bottom quark reduces more strongly than 
that of charm quark because of the heavier quark mass. 

From the radiation amplitude in Eqs.© and ([9]), the 
medium induced radiation probability to the first order 
in opacity can be expressed as 
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where C2, N, L and A± are, respectively, the Casmier in 
fundamental representation, the number, the thickness, 
and the transverse area of the targets, l g is the mean- 
free path of the gluon, a s = g 2 /4ir the strong coupling 
constant, P gq {x) = P(x)/x = [1 + (1 — x) 2 ]/x the splitting- 
function for q — > gq. |w(q_i_)| 2 is defined as the normalized 
distribution of momentum transfer from the scattering 
centers as shown in Ref. [T(|. The opacity factor L/l g 
reflecting the mean number of rescattering in the medium 
arises from the sum over the N distinct targets. 

The gluon formation factor inside Re(- • • ) in Eq. (fT8"|) 
reflects the destructive interference of the non-Abelian 
LPM effect [2(|. Comparing to the radiated gluon for- 



mation time t 
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j) for heavy quark jet 



goes through the static medium in Ref. [6[ , from the gluon 
formation factor in Eq. (fT8|) we see that the radiated gluon 
formation time in the medium with collective flow con- 
tracts by a factor 1 — v z because of the moving of the 
target parton. The pole shift of the quark propagator 



FIG. 2: The energy loss of charm, bottom and light quarks as 
function of flow velocity in the jet direction when E/fj, — 20, 
L/l g = 5. 

in the presence of collective flow leads to that the ra- 
diated gluon formation time in some interference pro- 
cesses changes by replacing uj m with to m — ljm as shown 
in the first term in Eq. (JTSJ) , so that the flow effects 
result in different gluon formation time for different in- 
terference processes. From Eq. (fT5)) we see clearly three 
different formation time, r/i = (1— v z )/{oj\+u> m ), 172 = 
(l—v z )l(L)i+u) m —u)M), Tf 3 = {1-v z )/uj m . The gluon 
formation factor should be averaged over the longitudi- 
nal target profile which is defined as (■ ■ • ) = J dzp(z) 
We take the target distribution as an exponential Gaus- 
sian form p(z) = exp(— z/L e )/L e with L e = L/2. 

As shown in Ref. Q , the scale of the hard scattering 
can be taken as Q 2 = 4E 2 , the kinematic limits of the 
gluon's transverse momentum are fi 2 < k\ < 4u>E. From 
Eq. (fT5|) we get the induced energy loss to the first order 
of the opacity as AE = j doJuidP 1 ^ 1 / 'dio . Its numerical 
results are shown in Fig. [2] for charm, bottom and light 
quarks as E / 'fj, = 20, L/l g = 5. At zero flow velocity, 
the light quark energy loss is larger than charm and bot- 
tom quark energy loss because of the dead cone from the 
mass effects of heavy quark jets as shown in Ref. @ 
in the static medium. However, the heavy quark energy 
loss increases with increasing the flow velocity because 
of the decrease of the dead cone by the flow effects in 
the medium with collective flow. The heavier the quark 
mass is, the more rapidly the heavy quark energy loss 
increases. The light quark energy loss decreases with in- 
creasing the flow velocity as shown in our previous work 
p^ |. In the region of 0.065<w z <0.13, the energy loss of 
light quark is a bit less than that of bottom quark but a 
bit larger than that of charm quark. In this region the 
average heavy and light quark energy loss is possible to 
be nearly the same. 

For A- A collisions at impact parameter b, with respect 
to collision number the average v z of rescattering centers 
in the QGP medium before freezeout can be expressed as 



where 



/ ^ Y I 27r / dr(v x cos(f) + v y sin(f)patAtB 
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) is the angle between jet and x axis, 



(19) 
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v y (\r+xiT — b/2|) are the flow ve- 



locity of the expanding elliptic medium along the minor 
and major semi-axes, n is the unit vector along jet direc- 
tion. Parton density p = p(r, b, r+nr), nuclear thickness 
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FIG. 3: The effective average energy loss of charm, bottom 
and light quarks as function of jet energy E in the presence of 
collective flow for — 10% central events of Au-Au collisions 
at RHIC energy. 

functions tA = tA(\r\), is = is(|r — b|). Cross section 
a = C a 2na 2 s / fi 2 [C a = 1 for qg and 9/4 for gg scatter- 
ing) obtained in 

pQCDdHH, P" = 9sT(t) at temperature 
T(r). Usin g th e data from a full 3D ideal hydrodynamic 
simulations [17j. we obtain (v z ) = 0.08 for 0— 10% central 
events of Au-Au collisions at RHIC energy, (v z ) is inside 
the region of 0.065<u z <0.13 as shown in Fig. [51 

The effective average energy loss of parton jet for A-A 
collisions can be written as 



JfrJ^Jd^Jd^ 
K ' Jd 2 rJdTt A t B 



^t A t B 



(20) 



The opacity L/l g in dP^ /d lo in Eq. (|18[) can be expressed 
as L/lg = J dTa(r)p(T, b, r + nr). Figj3] is the effective 
average energy loss of charm, bottom and light quarks as 
function of jet energy E in the presence of collective flow 
for 0—10% central events of Au-Au collisions at RHIC en- 
ergy. It is shown that the effective average energy loss of 
light quark is a bit less than that of bottom quark, but a 



bit larger than that of charm quark. The difference of the 
effective average energy loss among three quarks is very 
little, which implies that the light and heavy quarks have 
almost the same suppression of high px hadron spectrum 
for Au-Au collisions at RHIC energy. 

Conclusion — In summary, based on new potential in 
our previous work, we studied the flow effects on energy 
loss of the heavy quark jets in the quark-gluon medium 
with collective flow and give an explanation for the heavy 
quark energy loss puzzle. The collective flow changes the 
pole of the quark propagators, which lead to reduce sig- 
nificantly the dead cone of radiative gluon and increase 
obviously the energy loss of heavy quark jets. The col- 
lective flow also affect the gluon formation time because 
of the moving of the target parton and the shifting of the 
pole of the quark propagators. It has been shown that 
in the region of 0.065 < v z < 0.13, the light and heavy 
quark energy loss are nearly the same. For 0—10% central 
events of Au-Au collisions at RHIC energy, by using the 
data from 3D ideal hydrodynamic simulations we obtain 
that the average velocity of rescattering centers along the 
jet direction (i> 2 )=0.08, the difference of the effective av- 
erage energy loss among the charm, bottom and light 
quarks is very little. Our results shall have implications 
for comparisons between theory and experiment in the 
future. 
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